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Abstract 

We study the differential cross sections for electroweak gauge-boson and Higgs production at small and very small 
transverse-momentum qj . Large logarithms are resummed using soft-collinear effective theory. The collinear anomaly 
generates a non-perturbative scale q* , which protects the processes from receiving large long-distance hadronic con- 
tributions. A numerical comparison of our predictions with data on the transverse-momentum distribution in Z-boson 
production at the Tevatron and LHC is given. 



1. Drell-Yan-Like Processes 

Historically the Drell-Yan (DY) process [ 1 1 denoted 
the inclusive production of a virtual photon by quark- 
antiquark annihilation in hadron collisions and the sub- 
sequent decay into a lepton pair. Its main features are 
strongly coupled initial and color-neutral final states and 
so the photon case can easily be generalized to W- and 
Z-boson production. Even the Higgs production via 
gluon fusion can be described in a similar way. 
The transverse-momentum distribution of DY-like pro- 
cesses is one of the most basic observables at hadron 
colliders. It is used e.g. to extract the W-boson mass 
and width and is of great phenomenological relevance 
for Higgs -production at the LHC. Especially the regime 
of small transverse-momentum q 2 « M 2 is important, 
because it gives the largest contribution to the total cross 
section. Here qj denotes the transverse component of 
the boson 4-momentum q, while M 2 is its invariant mass 
q 2 . We thus consider: 

with q 2 = M 2 » q\ » A 2 QCD . 

dqr 

The hierarchy in this regime between the hard scale 
M and the collinear scale qj leads to large logarithms 
which spoil the perturbativity of fixed-order calcula- 
tions. These logarithms need to be resummed to all or- 
ders in perturbation theory to achieve a predictive result. 
Our approach O is to factorize the cross section using 
an effective field theory (EFT) and resum large loga- 
rithms via renormalization group (RG) techniques. The 
appropriate EFT to describe DY-like processes is the 
soft-collinear effective theory (SCET) [3 1, because it ac- 



counts for the complex structure of underlying scales 
originating from Sudakov double logarithms 0. 

2. Factorization using SCET 

SCET is an EFT of QCD. In general it describes 
any number of collinear modes, high energetic parti- 
cles (or Jets) with light-like momenta and soft modes, 
which mediate the only interactions between the differ- 
ent collinear fields. 
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Figure 1 : Momentum modes in SCET. 



In DY-like processes there are two collinear modes, 
defined by the two opposite light-like momenta of the 
colliding hadrons. The different momentum regions are 
best defined in lightcone coordinates. Therefore we in- 
troduce two light-like reference vectors n and h along 
the beam axis with n ■ h = 2. Now every 4-vector k 
can be decomposed into its collinear (k + ), anti-collinear 
(kJ) and perpendicular (k±) component, by projecting it 
onto n and h. 

The values of interest are the virtuality Vk 2 and the 
scalings of momenta: 



Scaling: k ~ (k + ,k_,k T ) with k\ 
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Figure 2: Scale uncertainties for different resuramation schemes. Figure 3: PDF uncertainties for different PDF sets. 



A virtuality of O(M) identifies the hard modes, which 
are integrated out like the produced DY-boson. The 
scaling is used to distinguish between the different 
collinear and soft modes (Figure [TJ. 

Up to power suppressed terms, the factorization using 
SCET leads to the following double differential cross 
section, where y denotes the rapidity of the DY-boson: 
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It consists of a hard function H, a sum over contributing 
partons and effective charges, a soft function W and two 
collinear functions S. The hard function contains the 
Wilson coefficients of the EFT. The soft function leads 
to scaleless integrals, thus does not contribute to all or- 
ders in perturbation theory: 



H(M,p) = \C(- 



-M 2 ,fi 2 



W = 1 +0(A 2 



Comparing the collinear functions S with the represen- 
tation of ordinary parton distribution functions (PDF) in 
SCET, it turns out they are just generalized xj depen- 
dent PDFs (gPDF): 
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Here the xj and p dependence is hidden in the logarithm 
L± = In (xj/i 2 ) . The Wilson coefficients are known, the 
soft corrections vanish and one can match the gPDFs on 
partonic level onto ordinary PDFs, only missing long- 



distance hadronic effects of ^(A^x 
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3. Collinear Anomaly and Resummation 

On the classical level the SCET Lagrangian re- 
spects the so-called rescaling symmetry. Since the 
two collinear fields can not interact with each other, 
each of their Lagrangians is invariant under the rescal- 
ing of momenta of the other one. At higher orders 
the collinear anomaly (CA) appears, the symmetry is 
broken by quantum corrections and restricted to joint 
rescaling, which introduces an unexpected invariant: 
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It turns out that this directly effects the matching of the 
gPDFs by generating a power-like dependence on the 
hard scale M: 

Si/nSj/n ^ (x 2 T M 2 ) F ^B i/N (£,L ± )B j/N (t,L ± ). 

This term ensures the RG invariance in the absence of 
soft contributions and is important for the resummation 
of large logarithms. 

The resummation of the hard function is simply done 
by using the RG equation: 



H(M,p) = H(M,p h ) ■ U(p h ,p) and set p h - 
Resumming the terms under the Fourier integral, 



M . 
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B i/N(£> L±)Bji N (g,L±) , 



is more subtle. It contains two types of logarithms, 
L± = In (xjp 2 ) from the collinear modes and In ^ T 
from the CA. Setting p ~ qj and expanding in a s leads 
to small L±, because xj is the conjugate variable of qj 
under the Fourier integral. At small qj this naive re- 
summation scheme leads to a large logarithm contained 
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Figure 4: Impact of hadronic effects on the intercept (top) and the 
peak-region (bottom) for Z-production. 
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Figure 5: Impact of hadronic effects on the intercept (top) and the 
peak-region (bottom) for Higgs-production. 



in 77 ~ a s In which spoils the perturbativity. 

To avoid this, our standard resummation scheme is to 

count 77 as 0( 1 ) and include higher order terms (in a s ) 

oftheCA. 

The standard resummation breaks down when 77 
reaches 1. This happens at the scale q*: 

q\ * 1.8 GeV, 7.7 GeV. 

To lower q T beyond q* one has to dismiss the demand 
of small L± by setting /j ~ q*. The appearing higher or- 
der terms of the CA form a Gaussian under the integral, 
which regulates it independently of qj, even at vanish- 
ing transverse-momentum q* » Aqqd > qr > 0. 

4. Uncertainties 

The first plot (Figure|2]l shows the renormalization 
scale uncertainties for Z-boson production at the Teva- 
tron. The error bands correspond to varying the default 
renormalization scale fid = It + <?f by a factor of two. 
The errors of the naive resummation (orange) lead to un- 
predictive results, because the terms of the CA are miss- 
ing. The standard resummation (red) gives considerable 
smaller error bands. The errors of the improved resum- 
mation (green) are somewhat smaller above q\, but sig- 
nificantly below. As a consequence all following plots 
are made using the improved resummation scheme. 

The second plot (Figure [3]) shows the PDF uncertain- 
ties for two different PDF sets at the Tevatron. The dif- 
ferent shape of the plot depends on showing ^- instead 



of ^p-. The first is used to point out the peak region, the 

latter for the intercept. The error bands correspond to 
one standard deviation to the center value. The uncer- 
tainties are around 5%, therefore lie within the renor- 
malization scale uncertainties. 

The hard function is independent of qj, thus can be 
regarded as an overall factor with constant uncertainties 
(Table [TJ. In the following plots the error bands corre- 
spond to the scale uncertainty. 

222 2 
4 = m z K = - m z 

"MX 1.000-;- I.334™ 
NNLL 1.0871*™ 1.131^ 
N 3 LL 1 H9+ 0006 1 130+ 0001 

Table 1: The hard function H (Mz,/j) at = Mz for space-like and 
time-like choices of . The uncertainties are obtained by varying ///, 
by a factor two about the default value. 

5. Long-Distance Hadronic Effects 

We model the non-perturbative effects with a Gaus- 
sian (blue) and a Dipole (red) factor in the gPDFs. The 
plots in Figure [4] show the impact of these effects on 
the intercept (top) and the peak region (bottom). By 
adjusting A NP we can fit the peak region onto experi- 
mental data without influencing the measurable rest of 
the cross section (the intercept can not be measured at 
hadron colliders). Since these effects should be univer- 
sal, we can set A^p in one measurement and use it as in- 
put for other distributions. All following plots are made 



3 




using the Gaussian model, because the differences be- 
tween the two models are marginal. 

Figure[5]shows the same plots in the Higgs case. The 
effects are strongly suppressed compared to the Z-boson 
production in Figure |4] The DY-like cross sections are 
protected from receiving large long-distance hadronic 
contributions by the CA. The effects scale with the ratio 
ANP /q*, which is much smaller in Higgs production. 



6. Final Results 



7. Conclusion 

As shown in the last plots, our approach of factoriz- 
ing the DY-like cross sections, using SCET and resum- 
ming large logarithms via RG-methods, leads to very 
good agreement of theory predictions and experimental 
data, together with small scale uncertainties. There have 
been a lot of approaches since the first resummation J6) 
in 1985, but this is the first time it was done directly in 
momentum space and it is free of Landau-pole singu- 
larities. Two important advantages of this approach are, 
it is straightforward to extend the calculation to higher 
orders in a s and A and the used methods are process in- 
dependent, therefore applicable to other problems [7 1. 



The last two plots show our final results, compar- 
isons with data on Z-boson transverse-momentum dis- 
tribution at CDF 0, Figure [6] without and Figure [7] 
with hadronic effects. Including these effects obviously 
improves the agreement of theory and data at small 
qr, while it does not influence the predictions above 
qr * 15 GeV. By using SCET we miss terms of 0(A 2 ), 
which become important at large q T . To receive a result 
for the whole g^-region, we match our result onto fixed- 
order calculations. The deviation at larger qj arises be- 
cause we only include matching at NLO fixed-order and 
should be reduced at NNLO. The matching correction is 
shown five times larger to make it visible. 
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